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Notwithstanding the visible maturity of the subject of Bose-Einstein Correlations (BEC), as 
witnessed nowadays, we would like to bring to ones attention two points, which apparently did not 
received attention they deserve: the problem of the choice of the form of C2(Q) correlation function 
when effects of partial coherence of the hadronizing source are to be included and the feasibility to 
model effects of Bose-Einstein statistics, in particular the BEC, by direct numerical simulations. 



I. INTRODUCTION 

The subject of Bose-Einstein Correlations 
(BEC) is so matured nowadays that (al- 
most) everything seems to be already answered 
and/or understood and what remains is just 
to systematically deduce from the experimental 
data information on the spatio-temporal struc- 
ture of the hadronizing source. That is the 
main reason of interest in BEC. Nevertheless, 
we would like to bring attention to two points 
which, in our opinion, are still worth of debate 
or, at least, worth to be remembered. They are: 

• What is the proper form of the two- 
particle correlation function C2{Q) in the 
case when one wants to account for the 
effects of the possible partial coherence of 
the hadronizing source 1]? 

• Is it possible to model numerically effects 
of Bose-Einstein statistics (BE), in par- 
ticular BEC, and in what way [2]? 

Since we address mainly readers already ac- 
quainted with the subject of BEC, no special 
introduction is offered. All necessary material 
can be found in [fl, I3| and references therein. 
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The above mentioned two points will be ad- 
dressed in two Sections that follow. We close 
with short summary. 

II. WHICH FORM OF CORRELATION 
FUNCTION IS THE CORRECT ONE? 

In most cases, when analyzing experimental 
data on BEC and discussing phenomenological 
models, one uses the following form for the two 
body correlation function {Q — \J — {pi — p2f 
and pi^2 are four-momenta of the observed iden- 
tical bosons): 

C2{Q)^l + \-n{Q-r). (Ill) 

Here is a 4-vector parameter, such that 
(r^)^ has dimension of length making the 
product Q ■ r ^ Qf^r^ ~ q dimensionless. It 
is usually regarded as representing dimension 
of the hadronizing source, but in reality it rep- 
resents the mean distance between the emission 
points of the two particles considered 0]. Pa- 
rameter A G (0,1), from the theoretical point 
of view, is usually understood as the degree of 
chaoticity of the source: totally chaotic source 
has A = 1 and totally coherent A = 0. Different 
shapes of the source were investigated in the 
literature: VL{Qr) = e"'?''; e-^''''; l/(l-hQr)2; 
[Ji{Qr)/{Qr)Y, to name a few [4[. The most 
advanced and complete discussion advocating 
this type of C2(Q) and justifying its structure 
can be found in [5|. 
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However, in a number of works it was 
strongly suggested that in the case when a 
hadronizing source is partially coherent the 
proper form of C2{Q) is the foUowing one [6]: 



C2(Q) = l + 2p{l~p)^n{Q-r)+p^n{Q-rf, 

(II.2) 

with r and defined as above, where p £ (0, 1) 
replaces A (retaining, however, essentially its 
meaning) . 

The imminent question arises then: which 
expression is the proper one? The answer was 
proposed in our paper [l| were we have shown 
that both expressions are correct in its own way, 
i.e., their form encodes information on the spe- 
cific features of hadronizing source which they 
describe; this point was not made apparent in 
the previous works 0, 0|- And so: 



Eq. describes situation in which 

hadronizing source can be regarded as 
consisting from the coherent and chaotic 
subsources acting independently in the 
proportion given by the chaoticity param- 
eter A 



• Eq. (jll.2[) describes situation in which 
there is only one hadronizing source but, 
for some reason, the phases of all particles 
are partially aligned. This can happen, 
for example, when hadronizing source is 
located in some constant external field, as 
was the case considered in 

To summarize: the choice of one or other form 
of C2(Q) presented here amounts to making 
a nontrivial assumption concerning the nature 
of the hadronizing source under investigation. 
This should be at least remembered and ac- 
knowledged, even if in practice only the first 
choice is nowadays used • 



III. HOW TO MODEL BEC 
NUMERICALLY 

The question of the numerical modelling of 
BEC is more important than usually antici- 
pated. The point is that to study multipar- 
ticle production processes one uses numerical 
simulations, the Monte Carlo event generators 
(MCEG) of different kinds. MCEG are based 



on classical probabilistic schemes whereas BEC 
is, by definition, quantum phenomenon and as 
such cannot be incorporated straightforwardly 
into a MCEG. The suggested cure was the use 
of the so called afterburners: one takes outcome 
of a given MCEG and changes accordingly mo- 
menta of the selected identical particles in such 
way as to fit the observed data on C2(Q). How- 
ever, it must be realized that by doing so one 
changes not only the original energy-momenta 
and/or multiplicities (for which one can correct 
later) but also (and usually unknowingly and 
in an unknown way) the physics of the model 
used as the basis of the MCEG chosen Q . 

The best solution would be to perform direct 
numerical simulations in which MCEG would 
start from the input containing already effects 
of BE statistics. What such input should be? 
The obvious suggestion is: the one possessing 
property that particles satisfying BE statistics 
tend to occupy in a maximal way the same 
state, i.e., they exhibit a bunching property. 
This property can be (at least in principle) 
modelled 0]. There exist already some exam- 
ples of such effort. In [10[ Metropolis method 



was used with fully symmetrized wave function 
to convert the set of N uniformly distributed 
identical particles into the set of N particles 
exhibiting the effect of BEC. Closer inspection 
shows that it happened because in this way 
particles were effecti vely bunched in the phase 
space forming what tlCI| called speckles. In 
one starts with single particle and, using re- 
jection method, adds to the iV-th particle the 
[N + l)-th one following the updated proba- 
bility as given by the fully symmetrized wave 
function for (TV -I- 1) particles. Again: the fi- 
nal distribution is characterized by bunches of 
identical particles in the phase space. In [l^ 
the main point was to account for the Negative 
Binomial (NB) character of the observed mul- 
tiplicity distributions P{N) by assuming that 
particles of the same charge are most likely be- 
ing produced in the same cells into which the 
phase space has been divided (it is rapidity in 
this work). All presented algorithms are very 
time consuming (especially the first two) and 
only the last one has been successfully applied 
to analysis of e+e^ data (and never again). 
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In Q we have summarized our effort aimed 
at improving this approach. From the exam- 
ples mentioned it is clear that the procedure 
of symmetrization of the initial set of identi- 
cal particles distributed somehow in the phase 
space takes too much time to be of any practical 
use. On the other hand, it leads to very inter- 
esting result, namely it shows that this proce- 
dure results in the effect of bunching of particles 
in some regions of phase space. Because bunch- 
ing is easier to simulate than symmetrization, 
it is this phenomenon which we propose to use 
as the cornerstone of the algorithm modelling 
BEG. Therefore we form bunches (called by us 
Elementary Emitting Cells - EECs) of particles 
in energy. It can be shown that using Bose- 
Einstein (or geometrical) form of distribution 
of particles in a single EEC one gets the charac- 
teristic proper BE form of < N{E) > together 
with the characteristic shape of the C2{Q) func- 
tion. When the original energy distribution is 
thermal-like (exponential with scale parameter 
T) then T is temperature seen in < N{E) >, 
whereas chemical potential present there is the 
main parameter describing BE distribution of 
particles in EEC The picture proposed re- 
sembles closely a quantum version of the clan 
model 13| (in which all particles in a clan are 



identical bosons). We call it therefore Quan- 
tum Clan Model (QCM) In this case one 
gets final multiplicity distribution in the form 
of Polya-Aeppli (geometric-Poisson) distribu- 



tion 
of 



14| , which differs from the NB distribution 



13| only for very small multiplicities. The 



strength of BEC, as given by parameter A in 
(jll.l[) . is very sensitive to the maximal allowed 
number of particles in the single EEC. There- 
fore, for extremely high multiplicity cases, A 
exceeds 2 (even for C2{Qj), a fact not noticed 
before. 



IV. SUMMARY AND CONCLUSIONS 



Let us summarize points raised here. 



by any external field which could make 
it partially coherent (depending on its 
strength). Therefore in situations where 
this cannot be assured its better to use eq. 
pi.2p . However, the trouble is that so far 
its form is elaborated only for two-particle 
BEC, multiparticle case still awaits its 
proper treatment. 

• The second point is that using one of 
the afterburners proposed in the litera- 
ture in order to change the outcome of 
the MCEG actually used, one accepts also 
(most times unknowingly!) all changes 
in the physical picture underlying this 
MCEG. The only way out would be to 
build a MCEG using the principle of 
BE statistics as its corner stone. It can 
be done by endowing MCEG from the 
very beginning of the numerical simu- 
lation process with property of bunch- 
ing of identical particles (via geometri- 
cal, or Bose-Einstein, particle distribu- 
tion in each bunch assured numerically). 
Effects of resonances, final state interac- 
tion and the like, can be (in principle) 
accounted for. They always reduce sig- 
nal of BEC. The most difficult problem 
one encounters is the implementation of 
corrections for nonconservation of energy- 
momenta and charge introduced during 
the Monte Carlo selection procedure used. 
We demonstrate that such program is 
possible but, at the moment, still far from 
the completion. 

Acknowledgments 

GW would like to acknowledge support ob- 
tained from the Bogolyubov-Infeld program in 
JINR and partial support of the Ministry of 
Science and Higher Education under contracts 
1P03B02230 and CERN/88/2006. 



• The first is that deciding, as it is usu- 
ally done, on the eq. (jll.ip when address- 
ing problem of BEC one tacitly assumes 
that hadronizing source is not infiuenced 
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